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We study the photon statistics properties of few-photon transport in an optomechanical system 
where an optomechanical cavity couples to two empty cavities. By analytically deriving the one- 
and two-photon currents in terms of a zero-time-delayed two-order correlation function, we show 
that a photon blockade can be achieved in both the single-photon strong-coupling regime and the 
single-photon weak-coupling regime due to the nonlinear interacting and multipath interference. 
Furthermore, our systems can be applied as a quantum optical diode, a single-photon source, and a 
quantum optical capacitor. It is shown that this the photon transport controlling devices based on 
photon antibunching does not require the stringent single-photon strong-coupling condition. Our 
results provide a promising platform for the coherent manipulation of optomechanics, which has 
potential applications for quantum information processing and quantum circuit realization. 
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I. INTRODUCTION 

The nonlinear effect is potential resource for quan¬ 
tum information processing [l|. For example, the photon 
blockade resulting from the nonlinearity is employed in 
single-photon (few-photonj transmission control |2| and 
optical state truncation |3[. Similarly, photon blockade 
is also an important feature in a lot of quantum device 
design such as fast two-qubit controlled-NOT gate 0, ef¬ 
ficient quantum repeaters ||, single-photon transistor Q 
and optical quantum computer [7j] . The rectifying device 
related with nonlinearity is the key device to informa¬ 
tion processing in integrated circuits [$]. Considerable 
efforts has been made to investigate the optical diodes 
0. Recently, various possible solid-state optical diodes 
have been proposed, for example the diodes from stan¬ 
dard bulk Faraday rotators [Tot integrated on a chip El, 
realized in opto-acoustic fiber |12| and from moving pho¬ 
tonic crystal El . A kind of optical diode based on pho¬ 
ton blockade effect also have been proposed, including 
photonic diode by a nonlinear-linear junction of coupled 
resonators 0 and optical diode of two semiconductor 
microcavities coupled via \^ nonlinearities (l5j. 

The nonlinear interaction between optical and me¬ 
chanical modes arising from radiation pressure force in 
optomechanical (OM) systems exhibit a lot of interest¬ 
ing nonlinear effects such as photon (phonon) blockade 
00 , optomechanical induced transparency [0 0 
and Kerr nonlinearity [0, (20] . Cavity optomechanics has 
receive d sig nificant attention both in fundamental exper¬ 
iments [0 10 ] and sensing applications 0 , 0 ] . Cur¬ 
rently, experimental technique of cavity optomechanics 
are still in the single-photon weak coupling regime [0 ] 
( g 2 < Kw m ), meanwhile it draws relatively few of works 
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as control devices in quantum information processing be¬ 
cause the prerequisite of the strong nonlinear is required 
0 . 0 ] . In order to utilize the nonlinearity of OM sys¬ 
tem in quantum information control, much attention has 
been paid to the photon blockade in OM system, in¬ 
cluding quadratically coupled OM systems 0, hybrid 
electro-optomechanical system [0 ], and ultrastrong op¬ 
tomechanics 1311. w here the strong coupling condition is 
required. Ref. [32]] has shown that strong photon anti¬ 
bunching can be achieved in two coupled cavities with 
weak Kerr nonlinearity, which motivate us try to achieve 
strong nonlinear effect in OM system in weak coupling 
regime. In this paper, we propose a scheme to realize an 
optical diode with optomechanical cavity coupled to two 
cavities. This scheme does not require the stringent con¬ 
dition that the single-photon optomechanical coupling 
strength g is on the order of the mechanical resonance 
frequency 0] or the coupling strength g is larger 
than the cavity decay rate k [291 ] . Our results show that 
photon blockade can be achieved both in strong and weak 
coupling regime because of the nonlinearity and the mul¬ 
tipath interference. By examining the second-order cor¬ 
relation function, rectifying factor 1Z and transport ef¬ 
ficiency T, we exhibit the characteristics of our system 
as photonic diode. Meanwhile, the single-photon trans¬ 
port can be controlled by the tuning of the frequency of 
the cavities in strong coupling regime. Surprisingly, by 
pumping two sides of the system (cavity L and R ), the 
device will embodies some characteristics like capacitor: 
photon storage-release (charge-discharge) and filtering of 
photon frequency. 


The paper is organized as follows. In Sec. II, we in¬ 
troduce the system and eigensystem of Hamiltonian. We 
discuss photons transport control in the cavity optome¬ 
chanical system as diode and capacitor in Sec. III. Dis¬ 
cussion and conclusions are given in Sec. IV. 
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FIG. 1: (Color online). Schematic of the cavity optomechan¬ 
ical system coupled with two cavities. Cavities L, R and C 
all can driven by the laser field with same frequency. 


II. MODEL AND HAMILTONIAN 

We consider a compound optomechanical system in 
which a cavity with a movable mirror is coupled with 
two cavities (L and R) with the coupling constant Jl and 
Jr, see Fig. [Q The system described by the Hamiltonian 
H = H sys + H purnp . By setting h = 1, the Hamiltonian 
H sys reads 


Hsys = E + Umb^b + g(tf + b)a} c ac 

j=L,C,R 

+(J L a^ L ac + Jua^ac + h.c.), (1) 

where ac, cll and cir are the annihilation operator for 
the photon mode of cavities G, L and R with frequency 
loc, wz, and lor, respectively, b is the phonon annihila¬ 
tion operator of the mechanical mode for the mirror with 
frequency cu m , g denote the coupling strength of radia¬ 
tion pressure. The cavity modes are driven by the laser 
with the same frequency wd, which can be described by 
H pU mp = J2j=L,c,R £ j( a ] e ~ ZU1Dt + hx - S >- In the rotatin S 
frame with Ho = c R^Ddjdj, we obtain 

Hs = E -vS + uj m b^b + g(b * + b)a' c ac 

j=L,C,R 

+(J L a' L ac + JrOr(ic + h.c.) 

+ E e f( a i + a j)> ( 2 ) 

j=L,C,R 

where A j = uij — u>d (j = L,C,R) are the detuning 
between the driving field and the jth cavity frequency, 
respectively. For this cascade configuration, cavity L and 
R are used as input and output ports in the side L and R. 
In this case, optomechanical cavity as an assisted-cavity, 
provides an intrinsically nonlinear interaction. 

We assume that the cavities (L and R) incoherently 
dissipate at rates n i (l = L, R) determined by the open¬ 
ness of the output channels and only classical driving 
fields are added to the quantum vacuum of the system, 
then according to the standard input-output relation [33|, 
the average output current (or photon stream) as num¬ 
ber of quanta emitted at time t from each cavity can be 
formally given by 

Qi(t) = m Tr[a\aip(t )], (l = L.R ), (3) 
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FIG. 2: (Color Online) The eigensystem of the Hamiltonian 
Hs in the zero-, one-, and two-photon cases, sub-area A, B 
and C denote multi-path interference in the system 


where p is the density operator of system. The evolution 
of the density operator p for the Hamiltonian Hr can be 
described by the master equation 

p = -i[H s ,p\+ E y Hla^p + ^(n th + l)V[b\p 

j=L,C,R 

+-^n th V[b^]p, ( 4 ) 

where nj and 7 are the cavity and mechanical energy 
decay rates, nth = [e-xp^m/ksTM — l )] -1 is the aver¬ 
age thermal occupancy number of the oscillator. T>[o] = 
2 oped — oi op — po^o is the Lindblad dissipation superop¬ 
erator. 

The eigen equation of the Hamiltonian H om = uj m b^b+ 
AcCL*a + g(tf + b)a^a can be expressed as 

Hom\s)c\n(s))m = E Sj „|s) c|™(s)} m, 
where the eigenvalues are 

E Sin = sA c + nu> m - s 2 S , 

2 

with S = , and the eigenstate 

|n(a)> = e 5 (fe - fct)/ "’"|n) (5) 

is the displaced number state. The eigensystem of the 
Hamiltonian Hr in the zero-, one-, and two-photon cases 
is shown in Fig. [2] We noticed that, the energy levels for 
optomechanical cavity (middle green line) will obtain a 
shift s 2 S caused by the nonlinear interacting with a fre¬ 
quency red (blue) detuning from the resonator resonance. 
This nonlinear shift can lead to bunched or antibunched 
photons in the OM cavity (the details are given in Sec. 
Ill C). This nonlinear effects also can appear in other 
cavities because of the coupling Jl and Jr. Especially 
in strong coupling regime g/n 1 , the system appears 
photon blockade, i.e. the probability for two photons in¬ 
side the cavity is largely suppressed due to the energy 
restriction. 
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The interference between multipath for two-photon ex¬ 
citation in cavities are partially responsible for the pho¬ 
ton antibunching effect shown in the sub-area of eigen- 
system diagram A, B and C. For area A, the two-photon 
in cavity L with state 12, 0,0)|n) m have two excitation 
path, one is direct excitation from low level in cavity L 
with state 11,0, 0)|n) m , the other is the tunnelling from 
OM cavity to left cavity with state |l,l,0)|n) m . The 
destructive interference between the two paths reduces 
the probability of two-photon excitation in the cavity. 
As well as area B and C. When the probability equal 
to zero, unconventional photon blockade [H, HHHU ap¬ 
pears in the cavity with no requirement to strong non¬ 
linear coupling coefficient g (even g/n < 1). Therefore, 
the compound optomechanical system can work as a sin¬ 
gle photon control device both in OM weak- and strong¬ 
coupling regime, which will be discussed in detail in next 
section. 


III. PHOTONS TRANSPORT CONTROL IN 
THE CAVITY OPTOMECHANICAL SYSTEM 


A. Optomechanical optical diode 


When the nonlinear effect for the right-going (k) is 
different from that for the left-going (— k) waves, i.e. the 
nonlinearity of the composite system is asymmetric, the 
rectification of one-dimensional photons transport can be 
controlled. The one way transport is called optical diode 
0. In this section, we will show that our compound 
system can worked as a photonic diode. 

We are interested in the statistic property of photons 
and its control. Usually, the frequency of the mechan¬ 
ical oscillator is larger than the strength of coupling of 
the radiation pressure, i.e., w m g. For simplicity, we 
can adiabatically eliminating the degree of the oscilla¬ 
tors. Including the decay rate of the cavities, we have 
non-Hermitian effective Hamiltonian as 


H eff = Y t( A i - iK j /2)a]a j + e^a] + aj)] 

j=L,C,R 

—i "Y^ -y-oja^ — Sa^ac — daQa^acac 

j=L,C,R 


+(J L a' L ac + JrOrClc + h.c.). 

We assume that the general state is 

IV4)> = C o {t)\0)+ Y Cj’W 0 ] I 0 ) 

j=L,C,R 

+ Y \ c n( t ) a \ a ] I 0 )- 

i,j=L,C,R Z 

Under weak pumping conditions, we have [32] 

C 0 » Cj > Cij. 


( 6 ) 


(7) 


Under this condition, one can obtain the steady state 
solution of the probability amplitudes, see Appendix. 


And the effect of quantum nonlinear features can be 
characterized by the second-order correlation function 
with zero-time delay. 




j = L,R,C. 


( 8 ) 


( 2 ) 

We notice that g - ^(0) < 1 indicates photon antibunching 

and g :j ; (0) > 1 indicates photon bunching, respectively. 
Antibunching corresponds to a reduced probability of two 
photons in the cavity at a given time, which is the oppo¬ 
site for bunching. The probability of two photons in the 
cavity will equal to zero if g ) ; (0) « 0 (photon blockade). 
For simplify, we set kl = Rr = Rc = k, ay = Ay — in /2 
(j = L,C,R). If a L = a R = a,e L = £ and £c,£r = 0, 
i.e., the system is only pumped on the left cavity with the 
magnitude e, the photon second-order correlation func¬ 
tion with no time-delay in left cavity can be obtained 


4 2) (o) 


4^2 + (Jl - JpKiF 2 

(.4 + J 2 r)k 2 - K,F 2 


X 


(j| + 4-Fi) 

(Jh-Fi) 2 


Fi\ 


(9) 


where K n = a+ac+n5 and F n = a(ac + nS) (n = 1,2). 
Setting A L = A R = Ac = A, we plot logarithmic g R \0) 
as a function of A and g in Fig. 3a. g R \o) ~ 0 repre¬ 
sents photon blockade, corresponding to the dark areas, 
which appears in two areas. The one is achieved in the 
low-right area in Fig. 3a with large values of coupling 
rate g , which means that the photon blockade is resulted 
from the nonlinear effects of radiation pressure. We call 
it conventional photon blockade (CPB). And the other 
appears in up-left with small values of g but strict lim¬ 
itations on other parameters, which means that it is re¬ 
sulted from the two-path interference. The interference 
between the two excitation path (the one exciton is from 
its own exciton, and the other is the jumping from its 
neighbor) is illustrated in Fig. [3jo. Because of the de¬ 
structive interference, the photon blockade phenomenon 
is also appearance, called unconventional photon block¬ 
ade (UPB). 

We now show the photon statistics properties and the 
controlling of the photon transport by comparing the an¬ 
alytical solution with the numerical results via solving 
the master equation 0. For the conventional photon 
blockade, the larger the ratio of g/n, the stronger the 
effect of blockade, shown in Fig. HJi. The correspond¬ 
ing detuning frequency can be derived from Eq. ([5]). As 
shown in Fig. [4}r, the strong photon antibunching can be 
obtained even if gjn < 1. 

In order to describe the characteristics of unidirectional 
energy transport. We define the rectifying factor 1Z and 
transport efficiency T as the normalized difference be¬ 
tween the output currents when the system is pumped 
through the left and right resonator (indicated by the 
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FIG. 3: (Color online) (a)The zero-time second-order corre¬ 
lation g^ (0) as a function of the coupling strength g and 
the driving detuning A (log scale). Other parameters are 
JL/uim = 0.5, Jr/lOtu = 0.01, n/cJm = 0.036. (b) The eigen- 
system of two excitation path interference in cavity L. 
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FIG. 4: (Color online) The equal-time second-order correla¬ 
tion g < l\ 0) as a function of the coupling strength g and the 
driving detuning A, set Al = A c = A, Ar = A + A lr- 
Other parameters are (a) = 0.5, Jr/uJui = 0.1, 

A LR = 0. (b) J L /uJm = 0.1, JR/u)m = 0.01, g / 0Jm = 0.01, 
k/ g = 1.3, A lr = 0.1 (other groups of parameters can be get 
by solving g^\ 0) = 0 in weak coupling regime). 


wave vectors k and — k, respectively) 


0 


iz = 


t l = 


r R = 


Q R [k) - QL[-k) 
Q R [k\ + Q L [-ky 
Q R [k] 

Qn[k] + Ql[/c] ’ 
Qd-k] 

Qr[— k) + Ql[~ k] 


( 10 ) 

( 11 ) 

( 12 ) 


7Z = — 1 indicates maximal rectification with enhanced 
transport to the left (left rectification), IZ = 0 indicates 
no rectification because = Ql [~&], while 1Z = 

+1 indicates maximal rectification with transport to the 
right (right rectification). In our system, cavity L and 
cavity R are both linear cavity. Therefore, there is no 
rectification (7 Z = 0) when only driving the left or right 



A/co 


FIG. 5: (Color online) The rectifying factor, transport effi¬ 
ciency and excitation number in cavity as a function of the 
driving detuning A, here we set Al = Ac = A, A_r = 
A -f- A LR- Other parameters are g/oj m — 5 x 10" 3 , Jl/oj m — 
5 x 10 -3 , Jr/ uJm = 5 x 10~ 3 , n/g = 0.2, A LR /uj m = 2 x 10” 2 . 


cavity (no asymmetric nonlinear effect). 

We discuss the rectification effect in conventional 
blockade regime (u> m > g > k ) shown in Fig. [ 5 ] When 
— is around —0.02, 1Z « + 1 which indicates that the 

UJm 

system allows photon transfer from left to right L —> R 
only, and the transfer is prohibited from right to left 
R —> L, the photon number from left-going (—k) field 
equal to zero in side L. Similarly, when i s around 
0.005, 7Z = — 1, which only allows the transport from 
right to left R —> L and NR^k) = 0. For is around 
-0.005, 1Z = 0, the photon number from left-going (— k) 
field equal to the photon number from right-going ( k ) 
field Nl{— k) = N R (k). If g < k < w m (unconven¬ 
tional blockade regime), shown in Fig. O one can obtain 
7Z ~ +1 when — = 0.022. We also can see 1Z = 0 for 

UJm 

■A. = 0.015, and TZ^-l for A_ = 0.003. Therefore we 
can conclude that no matter >g>K,oig<n< uj m 
by tuning the frequencies of the cavity R and L, one can 
adjust (or switch) rectification and two way transport. 


B. Single-photon source 

Photon blockade effect allows only single-photon trans¬ 
mission through the system. Now, we show that our 
device can work as single-photon sources. The system 
is only driven from left or right cavity, i.e., Sl = e, 
e R = Ec = 0 or er = £, el = £c = 0, the mean oc¬ 
cupation photon numbers 

N R (k) = N L (-k) (13) 

| JlJr£ 12 

JrO-l + &r{JI - a.L{ac + 6)) 

As shown in FigJTl the system only allowed single¬ 
photon transport no matter the light from left or right 
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FIG. 8: (Color Online) Pictorial representation and eigensys- 
tem of photon storage and release progress. 


FIG. 6: (Color online) The rectifying factor, transport effi¬ 
ciency and excitation number in cavity as a function of the 
driving detuning A, here we set A_l = Ac = A, Ah = 
A + Alr- Other parameters are g/uJm = 5 x 10“ 3 , n/g = 2, 
Jn/wm = 5x 10~ 3 , = 5x 10~ 3 , A LR /uj m = -2xl0 -2 . 
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FIG. 7: (Color online) The transport efficiency and second- 
order correlation of output ports as a function of the driving 
detuning A, here we set An = Ah = A, Ac = A + Anc- 
Other parameters are g/u} m = 0.2, Jn/w m = 0.1, Jr/uOtu = 
0.1, n/g = 0.1, Anc/wm = -1. 


when gjp(k) ~ g^\—k) ~ 0. The transport efficiency 
71 = 71 « 0.5, which means the output of system 
is single-photon state if the input is two-photon state. 
Under this condition, the device can control the single¬ 
photon transport in the channel or worked as a single¬ 
photon sources. This kind of device can only worked in 
strong coupling regime (g/n > 1) because there is no 
multi-path interference in output ports. We also notice 
that, 77 = 0 even if ^ cir,Jl ^ Jr- In order to 
achieve rectification, requires ec 0. 


C. Optomechanical optical capacitor 

As we have shown in Fig. [1] the photons in left cav¬ 
ity exhibit antibunching although the directly nonlinear 
interaction only appear in OM cavity. Similar to the non¬ 
linear shift 3d in OM cavity, the effective nonlinearity in 
cavity L and R can be equivalent to the resonance en¬ 
ergy shift Sl (Sr). If cavity L and R both appear photon 
antibunching effect due to the nonlinear shift and inter¬ 
ference while the OM cavity appears photon bunching 
effect, when we drive the cavity L and R , the photons 
can be stored in OM cavity. Reversing the process, one 
can release the photons. 

As shown in Fig. Ha) and (c), the system is a sym¬ 
metric structure. The field from left and right cavity can 
be regarded as input (+fc) of OM, meanwhile the field 
from OM cavity can be regarded as output (— k) of OM. 
In Fig. Hb), when {w L ,w R } > u D and {uj l ,ur} > uc, 
i.e. {A LC ,Arc} < 0, {A l ,Ah} > 0, the nonlinear fre¬ 
quency shift in left (right) cavity Sl(Sr) will enlarge the 
transition energy of two-photon excitation, which means 
the probability of two-photon state will be suppressed, 
the photon appears anti-bunching in cavity L (R ). At 
the same time, the nonlinear shift in OM cavity 3d will 
diminish the detuning between tunneling field uil (uir) 
and resonance frequency wc, the photon appears bunch¬ 
ing in OM cavity. Especially for ud ~ wc+3d, OM cavity 
exhibit strong bunching due to the resonance absorption. 
Under this condition, the probability amplitude of pho¬ 
tons in cavity OM will be much larger than in cavity L 
and R, photons can be stored in OM cavity. Reversing 
the process, as shown in Fig. Hd), when wc > ojd and 
loc > {wl,wr}, i.e. {A L c, Arc} > 0, A c > 0, the 
nonlinear frequency shift 3d will enlarge the transition 
energy of two-photon excitation, the photon appears an¬ 
tibunching in OM cavity. Meanwhile, the nonlinear shift 
in left and right cavity Sl and Sr will diminish the detun¬ 
ing between tunneling field ojc and resonance frequency 
the photon appears bunching in left and right 
cavity. Under this condition, the photons can be released 
from OM cavity. 
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FIG. 9: (Color online) The storage factor, storage-release 
efficiency and excitation number in cavity as a function of 
the driving detuning A, here we set A l = Ar = A, Ac = 
A + A lc- Other parameters are g/uJm = 5 x 10 —3 , n/g = 1, 
Jl/oOtu = Jr/ijOtu = lx IQ- 3 , ALc/oJm = —2 x 10' 2 . 


In order to describe the characteristics of energy stor¬ 
age and release. We define the storage factor 5 and 


storage-release efficiency M. 


5 = 

Qc[+k] — Ql[~ k] - Qr[- k] 

(14) 

Qci+k] + Ql[— k] + Qr\— fc] 

Ms = 

Qc[~\~k] 

(15) 

Qc[+^] + 

Mr = 

QL[-k] + QR[-k] 

(16) 

Qc[—k] + Ql[— fc] + Qr[~ k] 


where 5 = +1 indicates maximal storage with the en¬ 
hanced transport to the OM cavity, 5 = 0 indicates no 
storage and release, while 5 = — 1 indicates maximal re¬ 
lease with enhanced transport to the left and right cav¬ 
ities. Ms = 1 or Mr = 1 indicates the photons are 
totally stored in or released from OM cavity, respectively. 

For simplicity, assume that all parameters of cavity L 
and R are exactly the same in the following discussion. 
The photon number of the two cavity Qi[+fc] = Qn[+k] 
and Ql[~ k] = Qr[— k]. We discuss the storage effect 
shown in Fig. [9] When is around 0.02, S ss +1 
which indicates that system allows photon transfer into 
the OM cavity L — > C t— R only, and the transfer is 
prohibited out from the OM cavity L x— C — > R, the 
photon number in cavity L and R from OM cavity ap¬ 
proximately equal to zero. At this time, second-order 
correlation g y c {0 ) > 1, photons appear bunching effect 

in OM cavity, and (0) = 9l { 0) < 1, photons ap¬ 
pear antibunching effect in left and right cavity. The 
convergence filed (+fc) is bounded in OM cavity, system 
exhibits storage characteristic. As shown in Fig. EH sys¬ 
tem exhibits release characteristic. When — is around 

UJm 

0.007, S ~ — 1 which indicates that system allows pho¬ 
ton transfer out from OM cavity L C —> R only. The 


FIG. 10: (Color online) The storage factor, storage-release 
efficiency and excitation number in cavity as a function of 
the driving detuning A, here we set A l = Ar = A, Ac = 
A + A lc- Other parameters are g/uJm = 1 x 10 3 , k/ g = 2, 
Ji/Wi = J R / U)rn = 1 X 10 2 , ALc/ljJm = 2 X 1CF 2 . 


second-order correlation g^{0) < 1 , photons appear an¬ 
tibunching effect in OM cavity, and g^\ 0) = g^(0) > 1, 
photons appear bunching effect in left and right cavity. 
That is, divergent filed (—k) is released from OM cavity. 
We also notice that, when 5 = 1 and Ms = 1, indi¬ 
cates complete storage, no matter the field from left or 
right cavity can be stored in OM cavity, which is similar 
to capacitor charge process. While, when 5 = — 1 and 
Mr = 1 indicates complete release, the field in OM cav¬ 
ity can be released through the left and right cavity com¬ 
pletely, which is similar to capacitor discharge process. 
And the two progress can be controlled by the detuning 
of driving field A. On the other hand, like filter effect, 
there is no photon in the channel at the frequency which 
let 5 = 1 and Ms = 1 (complete absorption), but have 
no effect of the frequency which let 5 = —1 and Mr = 1 
(complete release). 

In the previous discussion, we ignore the effects of the 
mechanical thermal bath. Now, to investigate the influ¬ 
ence of the mechanical thermal temperature on the corre¬ 
lation function, we include the mechanical thermal reser¬ 
voir. Using master Eq.(4), in Fig. 11, we plot the mini¬ 
mum values of g\{ 0 ) as a function of the reservoir tem¬ 
perature, and the g\(0) versus affected by the ther¬ 
mal reservoir are also displayed in the inset. When the 
temperature below 1 mk (marked with the shadow area), 
the thermal heating nearly have no effects in Fig. EH a), 
because when the influence of the mechanical bath far be¬ 
low single-photon coupling rate, i.e. 7 nth g , the bath 
effect can be ignored. With current experimental tech¬ 
niques, one can easily set g/ 7 > 1 [36j, l37j, which means 
that a small value of phonon number nth can be toler¬ 
ance with no much effects. Also, we can clearly see that 
the antibunching effect becomes more and more weaker 
with temperature increasing. In UPB regime, as shown 
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FIG. 11: (Color online) The minimum equal-time second- 
order correlation g R \ 0 ) as a function of mechanical bath 
temperature Tm with mechanical frequency oj m = 0.1 GHz 
and dissipation rate 7 = lKHz. Other parameters are 
(a) A l = A r = Ac, g/oJm = 0.2, J L — J R = w m /10, 
K/bOm = 10 -2 . (b) A L = Ac, Alr/oJtu = 0 . 1 , g/uim — 10~ 2 , 
Jl/^tti = 0 . 1 , J R /uJm, = 0 . 01 , fc/g = 1.3. 

in Fig. fllT b). the antibunching effect is more sensitive to 
the bath temperature. And this quantum effect will dis¬ 
appear when the temperature over 5 mK (nth = 0.62). 
Fortunately, the current experiment conditions of ground 
state cooling can achieve n t h = 0.34±0.05 [22j. This pro¬ 
vides some ability to against the quantum decoherence of 
our system. Even so, to maintain the antibunching effect, 
the mechanical thermal noise still needs to be suppressed. 


IV. CONCLUSION 

In this paper, we employ the radiation pressure and 
the destructive interference effects to construct the con¬ 


troller of photon transport. By coupling the an cavity 
optomechanical system to two cavities, we show that the 
photon blockade can be achieved both in strong and weak 
coupling regime. In strong coupling regime, the photon 
blockade effects is mainly resulted from the nonlinear¬ 
ity of radiation pressure in optomechanical cavity, while 
in weak coupling regime, the photon blockade effects is 
mainly because of the interference between multipath for 
two-photon excitation in cavities. For few photon control 
of one-dimensional transmission, the system can worked 
as optical diode without the requirement of the strength 
of radiation pressure strong coupling, and the rectifica¬ 
tion of photons can be controlled by the detuning of 
driving field A. If we just drive the cavity from left 
or right cavity only, the system can function as single¬ 
photon source. Furthermore, when two fields transport 
into the OM cavity through the cavity L and R, the de¬ 
vice can store and release photons as a capacitor in ap¬ 
propriate parameter regime. These novel properties pro¬ 
vide a promising application of optomechanical system 
in quantum information processing and quantum circuit 
realization. 
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Appendix A: solution of the probability amplitudes 


We set kl = kr = kc = k, g, adiabatically eliminate the degree of the oscillators, we can obtain the equations 

of motion of the probability amplitudes under weak pumping regime Co C s i C s i S 2 - We drop higher-order terms 
in the zero- and one-photon probability amplitudes 


ICl = OLrCl-\- £lCq + JlCc, 

iCc — (otc + &)Cc + £cCo + JlCl + JrCr , 

iC‘R = Q!rCr + £rCq + JrCc, 

iClc = (oil + ac + 5 )Clc + £lCc + £cCl + V2Jl(Cll + Ccc) + JrClr, 

iClr = (oil + oir)Clr + £lCr + £rCl + JlCcr + JrClc, 

iCcr = (ocr + olc + S)Ccr + £cCr + £rCc + V2Jr(Crr + Ccc) + JlClr, 

iCll = 2ollCll + V%£lCl + V^JlClc, 

iCcc = (2 ac + 46)Ccc + 72ecCc + V2JlClc + V2JrCcr, 

iCRR = 2cxrCrr + V2£rCr + V2JrCcr , (Al) 


where = A l — in/ 2, ovr = A r — in/ 2, ac = A c — k/ 2 + <5. If we set the initial state is vacuum state, i.e. 
Co(0) = 1, C s i(0) = C' s i s2 (0) = 0,{sl,s2} e {L, C, R}, In the weak-driving regime, {£ C /Rc£l/kl,£r/Rr} < 1, the 
photon number is small, so we have Co(oo) « C'o(O), then the long-time solution of equations can be approximately 
obtained as, 


C L 

C c 

C R 

Cll 

Ccc 

Crr 


{—Jr + aR(ac + £>)]el + Jl(—olr£c + Jr£r) 

D[ ’ 

OtLOtR£c — JhOtR£L — JrOIlSR 

D[ ’ 

[ — Jr + Oic(otC + fi)]£R + Jr(—Oll£C + Jl^l) 

D[ ’ 

Cl 'l2j=(L,C,R) ^ L d £ j A JlCc Y^,j=(L,C,R) + JlJrCr Ylj=(L,C,R) 

7 W 2 ’ 

JlCl =(L,C,R) CL d £ j + Cc Ylj=(L,C,R) °CJ £ j + JrCr J2j=(L,C,R) C RJ £ j 

7Td 2 ’ 

JlJrCl J2j=(L,C,R) r L,j £ j + JrCc J2j=(L,C,R) r C,j £ j + Cr J2j=(L,C,R) r R,0 £ j 

7^7 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 


where the first term in Eq. (A5) describes two-photon state generated by driving field in cavity L , the second term 
describes two-photon excitation due to photon tunneling between OM cavity and left cavity with coupling rate Jl, 
the third term describes two-photon excitation due to photon tunneling between right cavity and left cavity through 









9 


the OM cavity, when the collective effect of this three progress let Crr ~ 0, photons exhibit blockade effect in cavity 
L. As well as Eq. (A6) and (A7). 


D i = JIcxl + JrOl r — a L a R (a c + 5), 

Di = 'y ] ot s [J 2 — a. s [a. s + cue + 5)][2J 2 {ot s + cue + 25) — a s (a; s + cue + 5){ac + 2(5)], 


s—L.R 


Il,l — Jl&r + [Jr — ( a L + otR)(ac + oil + (5)] [— a R (ae + ctR + 5 )(ckc + 2(5) + J R (ac + (a_jj + 2(5)] 

~ ctfi) — ot R (oie + a R (a L + ct#) + (3 q:r + cx R )5 + 2(5 2 + o;c(2aR + ctR + 3(5))], 

Il,c = d^R( a £ + ArK^c + + 25) — a fl (ac + 25)(—,/£ + (ay, + aj{)( a c + ^r + 5))], 

Il,r = Jr Jr[—J^(«c + or + 2(5) + + (ac + «R + <5)(ac + 2(5))], 

Ic,l = —J R {oiL + ctR){oLe + cur + 2(5) + afl(ac + 2(5) [—«/£ + (or + afl)( a c + ctr + <5))], 

lc,c = Jl[JrOil + JL a R — a R( a L + otR,)(ac + oir +6)\, 
lc,R = JlJr<X R (cXc + OIL + Or+ 2(5), 

Ir,l = Jr(&c + a R + 2(5) — c«r[</ 2 + (ac + cir + 6)(ac + 25)], 

Ir,c = JL&R(ac + oil + oir + 25), 

Ir,r = —JlJr(oic + &l + oir+26), 

cl,l = Jl[JrChl + J\olr — o.r(o.l + a R )(ac + oir + 5)], 

cl,c = ul[— JrCXl — J\olr + a R (aL + afl)(ac + «r + 5)], 

cl,r = Jr^ol — o:r(—J 2 + ai,(2ac + + ck/j + 25))], 

cc,l = Jlo:l[—JrO:l — J\olr + a R {ctL + a R ){ac + olr + 5)], 

ce,c = —Jt a R ~ a L[JR—(otL + atR)(ac + oll + 8)][Jr —oiR(ac + &R +$)] 

j cJ‘l[~ Jr(oil + oir) + c^R^r + cxe{ 2 aL + a^) + 0 ' R (a R + 5 ) + oil{ccr + 25 ))], 
cc,R = </rQ'r[—J^q:r — J 2 aR + 0 !R(aR + aR)(ac + or + 5 )], 
cr,r = Jl[JrOil ~ olr(—Jr + aR( 2 ac + «r + a R + 25 ))], 
cr.c = otR[—J R aL — JlC*r +atL(aL + atR)(ac + oil + S)\, 
cr,r = Jr[J^«r + — aR(aR + afl)(ac + ccr + 5 )], 

n,j = ^,j(^R ^ ^r,«r a R ),(i,j = {L, C,R}) (A8) 

The second order correlation functions with zero time-delay are 


si 2 ) (o) 


2|Crr| 2 


(|CR| 2 + |CR C | 2 + |CR fl | 2 +2|CRR| 2 ) 2 
2|Crr| 2 

|Cr | 4 ’ 



2|Ccc| 2 


(ic’d 2 + |Crc| 2 + \Ccr\ 2 + 2|Ccc| 2 ) 2 
2|C'c'c| 2 
|Cc| 4 ’ 



2 | C ' flfi | 2 


(|CA| 2 + |Ccr| 2 + |Crr| 2 + 2|C'_rh| 2 ) 2 
2|Crr| 2 
|Cr| 4 ■ 


(A9) 


The mean occupation numbers of three cavities are 


Nl = (|CR| 2 + |CRc| 2 + |CR fl | 2 + 2|CRR| 2 )iVo 

» |Cr| 2 5 V 0 , 


Nc = (|Cc| 2 + ICrcI 2 + ICcrI 2 + 2|Cc'c| 2 )-^ r o 

« |C c | 2 iVo, 
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N r = (\Cr\ 2 + \C C r\ 2 + \C lr \ 2 + 2 \C rr \ 2 )N 0 

« \C r \ 2 N 0 , (A10) 

where N 0 = (el/kl) 2 + (ec/«c) 2 + ( £r/kr ) 2 - r _ 



